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ABSTRACT 

We present an analysis of the diffuse ultraviolet (UV) background in a low latitude 
region near the Aquila Rift based on observations made by the Galaxy Evolution 
Explorer (GALEX). The UV background is at a level of about 2000 ph cm“^ s“^ sr“^ 
with no correlation with either the Galactic latitude or the 100 micron infrared 
(IR) emission. Rather, the UV emission falls off with distance from the bright B2 star 
HIP 88149, which is in the centre of the field. We have used a Monte Carlo model to 
derive an albedo of 0.6 - 0.7 in the UV with a phase function asymmetry factor (g) 
of 0.2 - 0.4. The value for the albedo is dependent on the dust distribution while g is 
determined by the extent of the halo. 
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1 INTRODUCTION 


The diffuse ultraviolet radiation field is comprised of sources 
ranging from airglow and zodiacal light (Murthy]|2014al to 


Lillie & Witt 


diffuse galactic light (Hayakawa et al. 1969 
|1976 ) with a possible extragalactic contribution at high 


Galactic latitudes (Henry 20021. Early observations in the 


field and their conclusions were reviewed by Bowyer (19911 


(20091. 


and Henry (19911 with a more recent review by Murthy 


There have been relatively few observations of the dif¬ 
fuse UV light at low latitudes, partly because many of the 
early observations were looking for an extragalactic back¬ 
ground (Sasseen et al. ]l9^ |Schiminovich et al.|20(iT| . One 
of the few surveys at low Galactic latitudes was the spectral 
imaging survey of the SPEAR {Spectroscopy of Plasma Evo¬ 
lution from Astrophysical Radiation: Edelstein et al. ( 2006|} 


satellite at wavelengths between 900 and 1750 A. 


Park et al. 


(20121 combined SPEAR and Galaxy Evolution Explorer 


[GALEX) data to study the diffuse UV radiation over a 30 
square degree region in the vicinity of the Aquila Rift. They 
deduced that most of the continuum emission was due to 
dust-scattered light but with contributions from molecular 
hydrogen fluorescence over part of the field. We have focused 
on a part of this field roughly corresponding to the region 
they call “Ophiuchus” and have used a Monte Garlo model 
of dust scattering to model the emission observed in the far 
ultraviolet (FUV) and the near ultraviolet (NUV). GALEX 
observed about 80% of the sky before it was turned off on 
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June 28, 2013 (Murthy 2014b I but had very few observations 


at low latitudes, primarily because of fears that the bright 
diffuse emission might damage the detectors. This constraint 
was relaxed late in the mission but the FUV power supply 
had failed by then. This region is the only location that in¬ 
cluded both FUV and NUV observations at a latitude below 
10 degrees. 


2 OBSERVATIONS AND DATA 

GALEX was a low Earth orbiting ultraviolet (UV) telescope 
under the NASA Small Explorer (SMEX) program whose 
primary goal was the observation of galaxies at low redshift. 
Along the way, it performed a survey of most of the sky 
in two UV passbands: the far ultraviolet (FUV: 1571 A) 
and the near ultraviolet (NUV: 2317 A). The spacecraft and 
the mission were described by [Martin et al.|(|2005[) and the 
calibration and data products by [Morrissey et al. (20071. 
The final data release (GR6/GR7) included all observations 
before the completion of the mission in 2013, barring a few 
proprietary observations taken as part of the CAUSE (Com¬ 


plete All-Sky UV Survey Extension) program. Bianchi et al. 


(20141 and Seibert et al. (20121 have combined all the in¬ 
dividual point sources detected in each observation into a 


single all-sky catalog and Murthy (2014b I has produced a 


diffuse map of the entire sky in each of the two bands. 

A GALEX observation consists of one or more visits 
pointed toward a single location in the sky. Because the vis¬ 
its may be on different days, each visit will have different 
contributions from the airglow and zodiacal light ( [Murthy] 
2014a I and so has to be processed independently if we want 
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Figure 1. FUV (left) and NUV (right) exposure time maps. The colour bars show the exposure times with black representing regions 
with no GALEX coverage. 



Figure 2. Original GALEX FUV image (left) and NUV image 
(right) for the lowest latitude observation, COROT_0507B. 


to extract the diffuse background. This was done by |Murti^ 
12014b I who ran all the GR6/GR7 visits through an au¬ 
tomated pipeline. He blanked out point sources using the 
merged catalog for each observation, binned the pixels into 
2' pixels and subtracted the foreground airglow and zodi¬ 
acal light using the formulation of Murthy (2014aI. The 
resultant data product consists of a separate data file for 
each visit containing only the astrophysical background and 
is available from the Mikulski Archive for Space Telescopes 

(mastQ 


This work comprises GALEX observations of an area 
of 15 square degrees centred at I = 30.0° and b — 10.0°. 
There were 227 observations with FUV data and 258 with 
NUV data, divided into 399 FUV visits and 567 NUV visits. 
Most of these observations were from the All-Sky Imaging 
Survey (AIS) with exposure times of 100 to 300 seconds but 
there were a few observations made by guest investigators 
with observation times of several thousand seconds each. 


^ https://archive.stsci.edu/prepds/uv-bkgd/ 


The exposure time per pixel of the co-added data is shown 
in Fig. For reference, we have shown the lowest latitude 
observation, which is also the deepest, at the full GALEX 
resolution in Fig. Note that there are very few UV stars 
in the field, even at these low latitudes. 

We binned the GALEX data in the field into 6' bins 
to increase the signal-to-noise in each of the two UV bands 
weighting each pixel by its exposure time. We have only 
used the central 0.5° radius of each observation to avoid edge 
effects in the data, as recommended by Murthy (2014b I. The 
resultant image in each of the two bands is shown in Fig. 
along with the lOOqm emission from the Infrared Astronomy 
Satellite (IRAS). There are two large exclusion zones in the 
UV data: in the centre around the star HIP 88149 (B2V, 
Vmag = 4.6, d = 207 pc); and at the lowest latitudes where 
the IR flux rises rapidly. Despite these gaps, it is apparent 
that there is no correlation between the UV and the IR 
intensities (Fig.|^. This is not surprising as the cross-section 


of the grains is much less in the IR than in the UV (Draine 


20031 implying that the thermal emission observed in the 


IR bands samples a much longer line of sight than the UV, 
which saturates over only a few hundred parsecs. 


3 MODELLING AND RESULTS 

We have developed a model for the dust scattered radiation 
over the entire sky | |Murthy|2015P and have used that model 
to predict the expected radiation in this region as a function 
of the Henyey-Greenstein ( [Henyey fc Greenstein|194T I opti¬ 
cal parameters: the albedo (a) and the phase function asym¬ 
metry factor (g). The albedo is the amount of light reflected 
by the grains while g (=< cos9 >) governs the angular de¬ 
pendence of the scattered radiation. Isotropic scattering is 
represented hy g — 0 with the amount of forward scattering 
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Figure 3. FUV (top left) and NUV (top right) diffuse maps of the 15°xl5° region centred at I = 30.0° and b = 10.0°. The IRAS 100pm 
map is at the bottom. There is no UV data in the centre of the region because of HIP 88149 (‘X’ in the map) and none at the lowest 
latitudes. 


increasing with increasing g. Although there are suggestions 


that the actual scattering is more complex (Draine 20031, 


the data have not been of sufficient quality to justify using 
more complicated scattering functions. 

The photons contributing to the diffuse background 
come from a relatively small number of O and B stars 


(Murthy & Henry 


1995 I for which a catalogue integration 


works well (|Henry 1977|. Here, we have used stars from the 


Hipparcos catalogue ( Perryman et al.|l997 l whose spectral 
type, distance, and B and V magnitudes are tabulated. The 
number of photons from each star at the desired wavelength 
were calculated using model spectra from|CasteIli fc Kurucz| 


| |2004[ ), as described by |Sujatha et al.| ( |2004[ ). Although stars 
from the entire Galaxy do contribute to the diffuse light, the 
diffuse flux in this field is dominated by a much smaller num¬ 
ber of stars (Fig. with 35% of the flux originating from 
just 10 stars and 65% from 100 stars. The brightest of these 
is HIP 88149 (66 Oph) which contributes 14% of the to¬ 
tal diffuse flux in the region and is responsible for the large 
exclusion zone in the centre of the UV observations. It is ap¬ 
parent from Fig.[^and Fig.|^that there is an enhancement of 
the diffuse flux in the UV — but no significant enhancement 
in the IR — around this exclusion zone. 

For the sake of this work, we have divided the sky into 
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Figure 4. GALEX FUV (left) and NUV (right) intensity plotted against IRAS lOOpm intensity. The y axis is in photon units (ph cm ^ 
s-l sr-l A-l). 
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Figure 5. The cumulative contribution to the flux from stars in 
the Galaxy shown. 


a 500x500x500 grid with a bin size of 2 pc on each side and 
filled each pixel with a hydrogen density corresponding to 
density of 1 cm“^ in the Galactic Plane with a scale height 
of 120 pc I Lockm^|1984 1 and a cavity of 50 pc in radius 
around the Sun ( [Welsh fc Lallement|2010| . We convert this 
into an optical depth for the dust using the “Milky Way” 


model of Draine (20031, which implicitly assu mes a uniform 
gas-to-dust ratio (N(H)/E(B-V)) of 5.8 x 10^^ (Bohlin et al. 


1978) throughout the Galaxy. More important to the scat¬ 


tered light in the observed region is the local distribution of 
dust as a function of distance. [Straizys et al.| ( |2003[ ) placed 
the bulk of the material in a cloud with its near edge at a 
distance of 225 ± 55 pc away from the Sun and a thickness 
of 80 pc. However, this would have placed most of the dust 
behind HIP 88149 which is at a distance of 207 pc from the 
Hipparcos catalogue ( Perryman et al.]|1997 ) with an E(B - 
V) of 0.19 magnitudes ( [Gudennavar et al.|2012p . We have 
therefore divided the dust into two clouds with a total red¬ 
dening in each pixel given by the all-sky map of |Schlegel| 
|et al.| ( |1998p . We placed the first cloud at a distance of 195 
pc from the Sun with a thickness of 10 pc and a maximum 


reddening of 0.2 magnitudes. The remaining dust, if any, was 
placed in the second cloud with a near edge of 205 pc and 
a thickness of 80 pc (Fig.j^. The spatial distribution of the 
modelled E(B - V) is in good agreement with the observed 
reddening (Fig. [^. 

We have run a series of simulations with different values 
of the optical constants to model the emission in the Aquila 
Rift. Because of time limitations, each run was for 50 million 
photons over the entire grid out of which only about 0.02% 
contributed to the flux in our field. At this stage in our 
modelling, it is difficult to match the data on a pixel by 
pixel scale and we have instead used the flux as a function 
of distance from HIP 88149 as a measure of the quality of 
the fit (Fig.[^. We have plotted the root mean square of the 
deviations of the model from the observed data as a function 


of albedo (a) in Fig. 10 and as a function of g in Fig. 11 for 


both the FUV and the NUV channels. The best fit model 
data implies an albedo (a) of 0.6 - 0.7 for both channels and 
a phase function asymmetry factor (g) of 0.2 - 0.4. 

Our current models are inefficient in that most of the 
computational time is spent in modelling the entire Galaxy 
and not just the Aquila Rift region. We have run much longer 
simulations (> 10® photons) for a more restricted set of op¬ 
tical constants and these are instructive in understanding 
the nature of the diffuse radiation. In each case, there is a 


large, bright halo around HIP 88149 (Fig. 131 the edges of 


which we see in the data. There is a trade-off in the mod¬ 
elled data between the albedo and g in that the level of the 
diffuse radiation increases with albedo, as expected, but de¬ 
creases with increasing g due to the sharpening of the halo. 
This is illustrated in the 3 ct contours for each of the FUV 
and the NUV in Fig. |12[ estimated using the procedure of 
[Lampton et al.| ( |1976[ ) . These limits do not take into account 
the uncertainty in the dust distribution. Although the bright 
scattering nebula around HIP 88149 is hidden in the exclu¬ 
sion zone of GALEX, because of the safety restrictions on 
observing near bright stars, it is prominent in the SPEAR 


data of Park et al. (2012) in the region they call “Ophi- 
uchus”. Note that there is no dependence of the modelled or 
the observed flux on the Galactic latitude. 

There are several uncertainties in our modelling, the 
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Figure 6. (a) GALEX FUV (b) GALEX NUV and (c) IRAS 
lOO^tm intensities plotted against angular distance from HIP 
88149. PU = ph cm-2 s"! sr”! A"!. 


Figure 9. Modelled FUV and NUV emission as a function of 
angular distance from HIP 88149 overplotted with the original 
FUV and NUV emission where the y-axis is in ph cm“^ s~^ sr~^ 
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Figure 7. Dust distribution as a function of distance to the cloud 
in a single line of sight. 
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Figure 10. Root mean square deviations as a function of albedo 
(a) for different g values. Lines connect points with identical g. 



Table 1. Previously determined a and g values 


a 


9 


Reference 


0.6 - 0.7 
0.28 ± 0.04 
0.40 ±0.10 
0.40 ±0.20 
0.45 ± 0.05 
0.45 ± 0.15 
0.50 

0.47- 0.70(136.2nm) 
0.55 - 0.72{176.9nm) 
0.42 
0.50 

0.60 ±0.05 


0.2 - 0.4 
0.61 ± 0.07 
0.55 ±0.25 

0.68 ±0.10 

0.90 


0.0-0.80 

0.44 

0.50 

0.75 ±0.15 
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Figure 8. Comparison of the modelled E(B - V) with the ob¬ 
served E(B - V) from Schlegel. 
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Figure 11. Root mean square deviations as a function of phase 
function asymmetry factor (g) for different a values. Lines connect 
points with identical a. 
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Figure 12. Estimated 3(7 confidence contours for a and g for 
FUV (solid line) and NUV (dashed line). 


most important of which is that we do not know the dis¬ 


tribution of the dust in any given line of sight (Seon et al. 


2011). Our modelling suggests that we require a substan¬ 


tial amount of dust near HIP 88149 in order to reproduce 
the observed halo around the star, yet we do not see a cor¬ 
responding infrared halo which would imply that the dust 
cannot be too near the star. As mentioned above, our cur¬ 
rent simulation wastes most of the photons in distant parts 
of the sky but because more than 10% of the diffuse flux in 



Figure 13. Modelled FUV image (left) and NUV image (right). 


the region is due to multiply scattered photons from outside 
the field of view, it is not trivial to reduce the scope of the 
simulation. This points to the immediate improvements to 
be made in our model. The first, as noted above, is to make 
the program more efficient by restricting the modelling to a 
much smaller area taking into account the fact that many 
of the scattered photons are non-local. The second is to im¬ 
prove the dust distribution. One important new input is the 
all-sky dust map based on PAN-STARRS and 2MASS pho¬ 
tometry by Green et al. (2015) which in combination with 
the limits put by the UV scattering observations may help to 
pin down the dust distribution and properties in this region. 


4 CONCLUSIONS 

We have used archival GALEX data combined with a Monte 
Carlo dust scattering model to probe the diffuse radiation 
in one field at low Galactic latitudes {b < 15°). Although 
the infrared 100 qm flux falls off rapidly from the Galac¬ 
tic Plane, there is no corresponding dependence of the UV 
flux on the latitude, at least in this restricted region near 
the Aquila Rift. Approximately 15% of the flux comes from 
HIP 88149, a B2 star in the field of view, which is respon¬ 
sible for an extended reflection nebula visible to a distance 
of more than 2° from the star. We fit the overall level of the 
diffuse background with an albedo of between 0.6 - 0.7 with 
a g between 0.2 - 0.4. The albedo is very dependent on the 
dust distribution in the line of sight and is therefore corre¬ 
spondingly uncertain. However, the value of g is constrained 
by the extent of the halo which would be too small if g were 
larger. These values are in reasonable agreement with other 
determinations of the albedo in a variety of regions (Table 
0 . 

Although the UV background follows a superficial cor¬ 
relation with the thermal IR emission over the entire Galaxy 


(Hamden et al.|2013 Murthy|2014b I 

there are clearly wide 

variations in the ratio, 
these were due to a dii 

Henry et al. 

(2015) suggested that 

Ferent, unknown contributor to the 


diffuse UV which was associated with the Galaxy but not 
part of the dust which contributes to the IR emission. They 
also pointed out that separate analyses of the background 
in small areas might individually be able to fit the flux by 
dust scattering but may not be consistent across regions. 

We find here that single stars may be the prime contrib¬ 
utors to the diffuse background over a wide area. Although 
at low Galactic latitudes, the scattered light is from a dust 
cloud at a distance of 200 pc from the Sun while the bulk 
of the IR emission is from an integrated column through 
much longer distances. Thus a complete understanding of 
the UV scattering depends on modelling both the diffuse 
emission over the entire sky as well as in individual loca¬ 
tions where local conditions are important. This has im¬ 
portant implications in the planning of observations with 
future UV missions such as the upcoming ASTROSAT mis¬ 
sion ( [Kumar et al.|2012[ ). The Galactic plane is excluded on 
safety grounds because of the presumed brightness of the 
UV background, as was done with GALEX. Our work here 
and with Voyager ( jMurthy et al.|1999| shows that there are 
dark regions even in the Galactic Plane. 
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